Absorbing State Phase Transition in presence of Conserved Continuous Local Field 
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We study absorbing state phase transition in one dimension in presence of a conserved continuous 
local field (CCLF) called energy. A pair of sites on a lattice is said to be active if one or both sites 
posses more energy than a pre-defined threshold. The active pair of sites are allowed to redistribute 
their energy following a stochastic rule. We show that, the CCLF model undergoes a continuous 
absorbing state transition when energy per site is decreased below a critical value. The critical 
exponents are found to be different from those of DP. 
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A state is called absorbing if it is impossible to leave 
the state. Existence of one or more absorbing configura- 
tions in a system raises a possibility of non-equilibrium 
phase transition between an active and an absorbing 
state Such absorbing state phase transitions (APT) 
are encountered in a variety of fields, which includes per- 
colation, spreading and chemical kinetics [3] ■ Numer- 
ous physical phenomena, like forest fire [Sl, epidemics 
[J, transport in random media synchronization Q 
and chromatography can be modeled as APT. The cor- 
responding critical behaviour forms a universality class of 
APT, formally known as the directed percolation (DP)- 
class, which has been realized Q convincingly in (2 + 1) 
dimension. 

Long ago, it has been conjectured by Grassberger and 
Janssen [8[ that in absence of any special symmetry or 
conservation law, a continuous phase transition into a 
single absorbing state governed by a fluctuating scalar 
order parameter, belongs to the DP universality class. 
Special symmetries, like particle- hole symmetry (9[, con- 
servation of parity [lO| . and symmetry between differ- 
ent absorbing states lead to different universalities. 
Also, conserved lattice gas model s II2II and conserved 
threshold transfer process (CTTP)[l3|, where the activ- 
ity held is coupled to the conserved density, show criti- 
cal behaviour different from DP. The DP-conjecture has 
raised the question, 'can systems with multiple absorbing 
states be in DP universality class?' In the pair contact 
process (PCP), initially introduced by Jensen [l^, the 
number of absorbing conflgurations grow to infinity with 
the system size. However, numerical simulations show 
that the critical behaviour is same as that of DP, which 
is further supported by phenomenological theories 15| . 
Other models, like the threshold transfer process la ], 
and dimer reaction models [T7| . which have infinitely 
many absorbing statcs(IMAS), also show static critical 
exponent same as DP, whereas the dynamical exponents 
(that characterize the spreading of localized perturba- 
tions) are non-universal [l^l; they depend on the nature 
of the absorbing states and initial conditions [l9| . Addi- 
tional conservation laws, like parity [20| . can change the 
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static exponents. 

More recent studies of Dickman and coworkers 
have renewed the interest in APT in systems having 
IMAS. They were able to show that the scaling behavior 
observed in sandpile models of self-organized criticality 
[2^ , is entirely governed by an underlying absorbing state 
transition into IMAS existing in equivalent fixed energy 
sand pile models (FES) [2^ . Corresponding universal be- 
haviour, which are different from DP, are attributed to 
the presence of coupling of the order parameter to the 
conserving height field [2J] . However certain specific per- 
turbations, like "stickiness" [25|, can drive these models 
to have critical behaviour same as DP. 

In all these models discussed above, the number of ab- 
sorbing states in a finite system is countable and grow ex- 
ponentially to infinity in the thermodynamic limit. Mod- 
els with continuous local field variable may show critical 
behaviour different from DP as (i) there is a possibility of 
having uncountably infinite number of absorbing states 
in these systems (even when th e sy stem size is finite) and 
(ii) as argued by Grassberger [2a], while for continuous 
variables one has "incomplete death" discrete systems do 
not allow such a process. However, the biological evolu- 
tion model (BEM) [23], which has a continuous dynam- 
ical variable, show DP-critical behaviour. Whereas the 
coupled map lattice models, where phase transition oc- 
curs to a synchronized (absorbing) phase, follow either 
DP or Karder-Parisi-Zhang universality class [2^, or a 
first order transition (29j depending on the non-linearity 
of the map [30| . Exactly which microscopic ingredients 
can make an absorbing state transition not belong to the 
DP class is an open and challenging problem. 

In this article we study a model in one dimensional 
lattice where the dynamical variables at the lattice sites, 
called energy, are continuous and their sum is conserved. 
They are updated pairwise, when one of the variable 
crosses a threshold value. We show that the system un- 
dergoes an absorbing state phase transition when the to- 
tal energy crosses a critical value. Study of several vari- 
ations of the models show that the critical behaviour of 
these systems is robust and forms a new universality class 
different from DP. 

The model is defined on a one dimensional lattice 
with periodic boundary condition with sites labeled by 
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i = 1,2...L. The dynamical variable called energy Ei, 
at each site i, is continuous and satisfies Ei = Ei^^ (pe- 
riodic boundary condition). A pair of neighbouring sites 
are called active when energy of one or both the sites 
cross a threshold value w. Otherwise, i.e. when both 
sites have energy less than w, they are inactive. Only 
the active pairs are allowed to exchange energy through 
an energy conserving dynamics, 
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where A is a parameter of the model and e is a random 
number distributed uniformly in (0, 1). Clearly the total 
energy E = ^^Ei is conserved by this dynamics. 

A special case of the model, w = corresponds to the 
kinetic model of markets introduced by Chakraborti and 
Chakrabarti (CC) [3l| in the literature of econophysics, 
where Ei is considered as wealth of an agent i who can 
interact with any other agent j (not necessarily its neigh- 
bour). Naturally, in this mean-field model the parame- 
ter A denotes the savings propensity. Variations of these 
models with variable savings propensity [3^ , provide the 
first explanation [33| 'why tails of the wealth distribution 
follow a power law called Pareto law [s^l'- Again, when 
A = 0, the distribution of energy (or wealth) follows a 
Gibbs distribution [s^ which has been observed in dis- 
tribution of income-tax return of individuals in several 
countries [s^. 

In fact A = case has been studied earlier [s^l on a 
lattice, in context of heat conduction, as a model of col- 
lisional dynamics of particles. First we consider this case 
A = with a finite threshold w (taken as w = 1, without 
loss of generality) . Note, that a different threshold w ^ 1 
would shift the transition point linearly. 

To study the properties of this model, particularly pos- 
sibility of phase transition and critical behaviour, we used 
Monte Carlo simulations. First let us define local activ- 
ity field Si = when the neighbouring sites i and i + 1 
are inactive, i.e. when both the sites have less than unit 
amount of energy. Or otherwise Si = \ (the bond join- 
ing i and i + 1 is active). Clearly, the activity (or the 
energy exchange) docs not die out when energy density 
e = Ei is larger than 1, as in this case the system 

has at least one active bond; starting from any arbitrary 
configuration, the density of active bonds p{t) = (s;), 
reaches a stationary value Ps as i cxd. Whereas for 
small energy density e <C 1, the activity is expected to 
die out, because it is highly improbable to have macro- 
scopic number of sites with energy larger than w. This 
indicates that there may be an absorbing state phase 
transition at ec < 1. 

In Fig. [1] we have plotted the steady state density of 
active sites ps as a function of energy density e. It is ev- 
ident that Ps takes a nonzero value as e crosses ec. Ps is 
plotted in log scale as a function of (e— ec) show a straight 
line by choosing the critical value Cc = 0.750(8) (see Fig. 
[IJb)). The orderparameter exponent /3 = 0.46(5), de- 
fined from the relation ps ^ [e — ec)^ , is clearly different 




FIG. 1. (Color online) The order parameter ps as a function of 
e. In log scale, ps is found to be linear in (e — Cc) by choosing 
ec = 0.750(8). Any small change (in steps of 5 = 0.0005) 
gives deviation from linearity (see inset (b)). Inset (a) shows 
that the orderparameter vanishes as ps ~ (e — ec)'', with 
/3 — 0.465, which is visibly different from Pcttp = 0.382. 
Here, the system size is L = 10*. 



(FiglUa)) from that of Pcttp = 0.382. 

One can obtain a few other exponents from the decay 
of p{t) from a fully active configuratior0 with p(0) = 1. 
After an initial decay p{t) ^ it approaches the steady 
state value ps in the t ^ oo limit. So p must scale as 



p(f,e) = t-"^(ii/''ii(e-ec)) 



(2) 



Thus, one expects that p{t) for different values of e 
(shown in the Fig. ^ collapse into a single scaling func- 
tion T, when t"p{t) is plotted against t\e — ed^n in log 
scale. The main figure here shows the data collapse when 
we use a = 0.19(5), = 2.64. Since at the critical point 
p{t, Cc) ~ t~"J^(0), one can obtain both Cc and a directly. 
Resulting Cc and a are consistent with those obtained 
from the data collapse. Again, in t — cxd limit, p vanishes 
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This can happen only when J-{x) 



a = /3/j/||. 



Since all three exponents /3, a and are calculated in- 
dependently, one can check if the above scaling relation 
holds. In fact, it holds to a great accuracy for the values 
of /?, q; and calculated here. 

Now we turn our attention to the finite size scaling of 
p{t) at the critical point. Again, from a fully active state 
p{t, L) decays as t~", indicating a scaling form 



p{t,L)^t--g{tiL^), 



(3) 



where z is the dynamic exponent. Thus, one expects that 
p{t) for different values of L collapse to a single function 



^ This can be generated by taking an initial condition such that 
all sites of one sub-lattice have energy larger than unity. 
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FIG. 2. (Color online) Inset : Decay of p{t) from a fully active 
state to the steady state value pa is shown in the inset for L — 
1000 and e = 0.65,0.67,0.69,0.75,0.77,0.8, and 0.85 . Main 
figure : Data collapse for off-critical simulations according to 
the scaling form Eq. ([2]) is obtained for a = 0.19(5) and 
i^ll = 2.64. 











L=128 

L=200 


1 








L=350 

V L=512 












_ 0.01 

p 










; 0.0001 


1 1 . .i.i.l 








100 f 10000 




le+06 
1 .... 





FIG. 3. (Color online) The decay of p{t) at = 0.750 from 
a fully active state, averaged over 10* ensembles, is shown in 
the inset for different L = 128, 200, 256, 350, and 512. A line 
with slope a = 0.195 is drawn to emphasize that p{t) ~ t~" . 
The main figure shows the data collapse for finite size scaling 
Eq. @ using z = 1.381. 



when plotted against t/L^. This is described in Fig. O 
The inset there shows variation of p{t) for different L = 
128, 200, 256, 350, 512, which were made to collapse to a 
single function using a = 0.19(5) and z ~ 1.38(1). Now, 
assuming the scaling relation z = ly^^/v^, one can obtain 
v± = fj/az. 

The exponents we obtained using Monte-Carlo simu- 
lations are listed in Table HI along with the exponents of 
DP, BEM and CTTP. Clearly, the exponents of CCLF 
model are quite different from the DP. Even the static 
exponents, which are believed [l^ to be same as that of 
DP for models with multiple absorbing sates, are now 
different. Though some of the exponents are close to the 
same obtained for CTTP, particularly the order param- 
eter exponent /3, and the spreading exponent a are very 
different. Since the value of /3 crucially relies on the esti- 
mation of the critical point, we provide a careful study of 
Cc based on Ref. [1^ in Appendix-I. From this analysis it 
is concluded, beyond reasonable doubts, that /3 = 0.46(5) 
and thus the phase transition in CCLF model belong to 
a new universality class. 

One can possibly reason it to the existence of a con- 
tinuous conserved field. Further study in this direction 
is required to identify, what made this absorbing state 
phase transition different from the usual ones. 

TABLE I. Critical exponents of DP, BEM and CCLF model. 

P v\\ a z 

DP [2] 0.276 1.733 0.159 1.581 

BEM [27] 0.276 1.067 .259 1.364 

CCLF 0.46(5) 2.6(4) 0.19(5) 1.38(1) 

CTTP [13] 0.382 2.45 0.141 1.393 



In the following we discuss some possible directions 
of studies which may explain the new universality. The 
model studied here, is quite similar to the sandpile mod- 



els [22| of self organized criticality, except that the local 
field variables in CCLF model are continuous. It has 
been pointed out in Ref. [2^, that although the critical 
behaviour of sandpile models crucially depends on the 
details of the dynamics and the spatial dimension, they 
are unstable to certain specific perturbation (called stick- 
iness) and flow generically to the DP universality class. 
It would be interesting to ask, if the model studied here 
is unstable to perturbations. Some generic variations, 
which may change the critical behaviour of CCLF, are 
discussed below. 

The bulk dynamics of CC model (A ^ 0) does not 
satisfy detailed balance, as for any arbitrary (Ei,Ej) — 
{El,E'j) it is not possible to have {E-,E'j) {Ei,Ej). 
Thus, A ^ is a singular perturbation and it may change 
the universality class. The quenched disorder, introduced 
by taking a distributed savings propensity {A^} [3^ may 
also change the critical behaviour. 

In conclusion, we have studied absorbing state phase 
transition in a model with conserved continuous local 
field (CCLF). In one dimension, the model is defined on 
a lattice with sites i having continuous variable Ei called 
energy. Two neighbouring sites i and i + I can exchange 
energy when one or both sites have energy larger than a 
predefined threshold w. The exchan ge d ynamics is sim- 
ilar to the wealth exchange model [31| studied earlier, 
by considering Ei as wealth of the agent i. Clearly, the 
system is active {i. e. surely some of the neighbour- 
ing sites keep exchanging their wealth) when the energy 
density e is larger than w (set to be unity, without loss 
of generality). Whereas for e ^ 1 it is improbable to 
get macroscopic number of sites which have wealth larger 
than w indicating that the system falls into one of the un- 
countably infinite number of absorbing states. We show 
that, CCLF model undergoes a continuous phase transi- 
tion at Cc = 0.750(8) with critical exponents (3 = 0.46(5), 
i^ll = 2.6(4), a = 0.19(5), and z = 1.38(1), which are 
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very different from those of DP. However earlier stud- 
ies [isj . botli numerical and analytical, have shown that 
the static exponents of absorbing state transition to in- 
finitely many absorbing states are same as those of DP, 
whereas dynamic exponents differ. It is surprising, that 
in CCLF model, even the static exponents differ from 
DP and form a new universality class. One may argue 
that the existence of the conserved field is responsible for 
this behaviour. But, this can not be the sole argument 
as recent studies [2^ show that it is possible to get DP 
behaviour in presence of conserved fLeld(s). 

Appendix-I 

The numerical estimation of the critical exponents of 
CCLF model are listed in Table-I. The critical exponents 
i^ll and z are quite close to that of CTTP, whereas /3 and a 
are substantially different. The main claim of this article, 
that CCLF model form a new universality class, is based 
on this difference. Since the value of /3 strongly depends 
on the the critical point , here we estimate Cc through a 
careful and systematic numerical analysis based on Rcf. 
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FIG. 4. (Color online) (a) p versus t for e = 0.76 and 0.765 are 
fitted with the scaling form ((4|. (b) pt"" — 1 as a function 
t in log scale is linear for both e = 0.76 and 0.765 when A = 
2.64; the solid lines indicates a power- law i^/^ll with v\\ = 2.64 
and coefficients 0.0178 and 0.0269 respectively. From these 
coefiicients we find that B = 1.82 and = 0.75022. The 
inset of (a) shows that p{t,ec) obtained from p{t,e — 0.76) 
using Eq. (O varies as t~°' with a = 0.19, when Cc is taken 
as 0.7503. Here L = 1024. 

Near the critical point, the order parameter p{t, e) fol- 
low Eq. ©; thus 



p{t, e) = + B{e ~ Cc) i^/"^" ) 

~p(t,e,)(l + B(e-e,) t^l^w 



(4) 



where we have used the Taylor's expansion of scaling 
function T(x) at a; = upto the first term, i.e. J'{x) = 
A{1 + Bx). Clearly, this functional form is valid for 
t < |e-e,|"''ii. 

Using this functional form, which is valid for < <C |e — 
ecl^"^", one can obtain by suitably choosing A such 
that e)t" — 1, as a function of t, is linear in log 

scale; the slope and the y-intercept determine the v\\ and 
B{e — Cc) respectively. For a system of size L = 1024, 
p is shown in Fig. |4] (a) for two different values of e = 
0.76 and 0.765. Figure IDJb) shows that A~^p{t,e)t'^ - 



1 is linear in log scale for both values of e, with a = 
0.19 and A-^ = 2.64. The slope of both the straight- 
lines turns out to be l/;/|| = 0.379, which is consistent 
with j^ii = 2.6(4) obtained earlier from the data-collapse. 
From their y-intcrcepts wc get, B{e — ec) = 0.0178, 0.0269 
respectively for e = 0.76 and 0.765, which provide Cc = 
0.75022 and B = 1.82. Using these parameters, we have 
calculated p{t, e) as a function of < for e = 0.76, 0.765 
(shown as solid lines in Fig. IH^a)). The discrepancy for 
large t, in case of e = 0.765 is due to the fact that we have 
approximated J-{x) in Eq. ^ up to the linear order. 

A better estimation of Cc can be done now by using 
Eq. (HD, 



p{t,e) 



\ + B{e~-ec) ^^/''ii 



(5) 



Since we know the value of B and v\\ , Cc can be used as a 
fitting parameter in the above equation such that p{t, Cc) 
obtained from the numerical values of p(t, e) using Eq. 
^ shows a power-law. The inset of Fig. S] (a) shows 
that p{t,ec), obtained from p{t,e = 0.76) by choosing 
Ec — 0.7503, decays algebraically for four decades with 
exponent a = 0.19. Thus the final estimate of the critical 
point for L ^ 1024 is = 0.7503. 
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FIG. 5. (Color online) (a) ps as a function of e for different 
L = 128,256,512,1024 and lO". Difference in p^ is visible 
only after sufiicient magnification, (b) pa for different L are 
merged onto the same for L = 10* by using Eq. ((6|. The 
inset here show variation of ec{L) with 1/L, which is linear. 
The critical value obtained for L = 1024 following Ref [i^l is 
shown as a circle. The asymptotic value of the critical point, 
extrapolated for L oo is Cc = 0.7508. 

In general, the order parameter and thus the criti- 
cal point depends strongly on system size L. In Fig. 
[5] (a), wc have plotted ps as a function of e, for L = 
128, 256, 512, 1024 and 10''. These curves do not show 
any significant dependence on L. The small difference 
can be adjusted by shifting the x-axis by a small amount 
6l, which results in 



Psie) 



pTie + SL)-{e-e,iL)f, 



(6) 



where ec{L) = ec(oo) — 5^. Wc calculate edL) as the 
value, where the log scale plot of ps versus e — ec{L) is 
linear. In fact, with these choice of ec{L), all the different 
curves of p^ could be merged on to the curve for L = 10^ 
(see Fig. E] (b)). The inset of Fig. M (b) shows the 



5 



small variation of ec{L) with which asymptotically 
approach to Cc = 0.7508 as L cxd. Again, the value 
of Cc obtained in this Appendix for L = 1024, using a 
method described in Rcf. consistently fall on this 

curve (denoted by a circle in the inset of Fig. E] (b)) . 



These analysis suggest that the estimates of Cc and 
the critical exponents /3 and a do not change apprecia- 
bly with system size. The current values of exponents 
obtained for system size L = 10"* are sufficient to con- 
clude that the critical behaviour in CCLF model belongs 
a universality class different from CTTP. 



[1] J. Marro and R. Dickman, Non-equilibrium phase tran- 
sitions in lattice models (Cambridge University Press, 
Cambridge, 1999); 

[2] H. Hinrichsen, Non-equilibrium Critical Phenomena and 
Phase Transitions into Absorbing States, Advances in 
Physics 49, 815 (2000); M. Henkel, H. Hinrichsen, and 
S. Liibeck, Non-Equilibrium Phase Transitions, Berlin: 
Springer (2008). 

[3] E. V. Albano, J. Phys. A 27, L881 (1994). 

[4] D. MoUison, J. Roy. Stat. Soc. B 39, 283 (1977). 

[5] S. Havlin and D. ben Avraham, Adv. Phys. 36, 695 
(1987); J. P. Bouchad and A. Georges, Phys. Rep. 195, 
127 (1990). 

[6] P. Grassberger, Phys. Rev. E59, R2520 (1999); V. Ahlers 
and A. Pikovsky, Phys. Rev. Lett. 88, 254101 (2002). 

[7] K. A. Takeuchi, M. Kuroda, H. Chate, and M. Sano Phys. 
Rev. E 80, 051116 (2009); H. Hinrichsen, Physics 2, 96 
(2009). 

[8] P. Grassberger, Z. Phys. B 47, 365 (1982); H. K. Janssen, 

Z. Phys. B 42, 151 (1981). 
[9] J. W. Essam, J. Phys. A 22, 4927 (1989). 
[10] I. Jensen, J. Phys. A 26, 3921 (1993). 
[11] H. Hinrichsen, Phys. Rev. E 55, 219 (1997). 
[12] M. Rossi, R. Pastor-Satorras, and A. Vespignani, Phys. 

Rev. Lett. 85, 1803 (2000). 
[13] S. Lubeck and P. C. Heger, Phys. Rev. E 68, 56102 
(2003). 

[14] I. Jensen, Phys. Rev. Lett. 70, 1465, (1993); L Jensen, 

Phys. Rev. E 47, Rl (1993); L Jensen and R. Dickman, 

Phys. Rev. E 48, 1710 (1993). 
[15] M. A. Munoz, G. Grinstein, R. Dickman and R. Livi, 

Phys. Rev. Lett. 76, 451,(1996). 
[16] J. F. F. Mendes, R. Dickman, M. Henkel, and M. C. 

Marques, J. Phys. A 27, 3019 (1994). 

[17] J. Kohler and D. ben-Avraham, J. Phys. A 24, L621 
(1991); D. ben- Avraham and J. Kohler, J. Stat. Phys. 
65, 839 (1992); E. V. Albano, J. Phys. A 25, 2557 (1992); 
A. Maltz and E. V. Albano, Surf. Sci. 277, 414 (1992). 

[18] M. A. Munoz, G. Grinstein, R. Dickman and R. Livi, 
Physica D 103, 485 (1997); C. Lopez and M. A. Munoz, 
Phys. Rev. E 56, 4864, (1997). M. A. Munoz, G. Grin- 
stein and R. Dickman, J. Stat. Phys. 91, 541 (1998); F. 
van Wijland, Braz. J. Phys. 33, 551 (2003). 

[19] G. Odor, J. F. Mendes, M. A. Santos, and M. C. Marques, 
Phys. Rev. E 58, 7020 (1998). 

[20] M. C. Marques and J. F. F. Mendes, Eur. Phys. J. B 
12, 123 (1999); J.F.F. Mendes, Braz. J. Phys. 30, 105 
(2000). 

[21] R. Dickman, A. Vespignani, and S. Zapperi, Phys. Rev. 



E 57, 5095 (1998); A. Vespignani, R. Dickman, M. A. 
Munoz, and S. Zapperi, Phys. Rev. Lett. 81, 5676 (1998). 

[22] H. J. Jensen, Self-Organized Criticality (Cambridge Uni- 
versity Press, Cambridge, 1998); D. Dhar, Physica A 
369, 29 (2006); A. Ben-Hur and O. Biham, Phys. Rev. 
E 53, R1317 (1996). 

[23] This relation has recently been put into doubt, see, P. K. 
Moha nty and D. Dhar, Physica A 384, 34(2007); S.-C. 
Park, \arXiv:1001.33M 

[24] A. Vespignani, R. Dickman, M. A. Munoz, and S. Zap- 
peri, Phys. Rev. Lett. 81, 5676 (1998); M. A. Munoz, R. 
Dickman, A. Vespignani and S. Zapperi, Phys. Rev. E 
59, 6175 (1999); R. Dickman, M. A. Munoz, A. Vespig- 
nani and S. Zaperi, Braz. J. Phys. 30 27 (2000). 

[25] P. K. Mohanty, Deepak Dhar, Phys. Rev. Lett. 89, 
104303 (2002). 

[26] P. Grassberger, Phys. Rev. E 59, R2520 (1999); P. Grass- 
berger, In Nonlinearities in Complex Systems, edited by 
S. Puri and S. Dattagupta (Narosa Publishing House, 
New Delhi, 1997), p. 61. 

[27] A. Lipowski and M. Lopata, Phys. Rev. E 60, 1516 
(1999). 

[28] M. Kardar, G. Parisi, and Y. C. Zhang , Phys. Rev. Lett. 

56, 889 (1986). 
[29] P. K. Mohanty, Phys. Rev. E 70, 045202(R) (2004) 
[30] A. S. Pikovsky and J. Kurths, Phys. Rev. E 49, 898 

(1994); L. Baroni, R. Livi, and A. Torcini, Phys. Rev. E 

63, 036226 (2001); M. A. Munoz and R. Pastor-Satorras, 

Phys. Rev. Lett. 90, 204101 (2003). 
[31] A. Chakraborti and B. K. Chakrabarti, Eur. Phys. J. B 

17, 167 (2000). 

[32] A. Chatterjee, B. K. Chakrabarti and S.S. Manna, Phys- 
ica Scripta T 106, 36 (2003); ihid Physica A 335, 
155(2004). 

[33] P. K. Mohanty, Phys. Rev. E 74, 011117 (2006); Urna 
Basu and P. K. Mohanty, Eur. Phys. J. B 65, 585 (2008). 

[34] Pareto V., Le Cours dEconomie Politique (Macmillan, 
London) 1897. 

[35] A. A. Dragulescu, and V. M. Yakovenko, Eur. Phys. J. 

B 17, 723 (2000). 
[36] A. C. Silva and V. M. Yakovenko, Europhys. Lett. 69, 

304 (2005). 

[37] C. Kipnis, C. Marchioro, and E. Presutti, J. Stat. Phys. 
27, 65 (1982). 

[38] P. Bak, C.Tang and K. Wiesenfeld, Phys. Rev. Lett. 59, 

381 (1987); J. Phys. A 38, 364 (1988); D. Dhar and R. 

Ramaswamy, Phys. Rev. Lett. 63, 1659 (1989). 
[39] S. S. Manna, Physica A 179, 249(1991); S. S. Manna, J. 

Phys. A 24, L363 (1991). 
[40] E. Perlsman and S. Havlin, Europhys. Lett. 58, 176 

(2002). 



